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1 . INTRODUCTION 


Linear wave propagation in hot collisionless plasmas is described 
by the linearized Vlasov and Maxwell equations. In uniform media, the 
utilization of spatial and temporal transforms of those equations leads 
to the consideration of integrals of the type (Hilbert transform) 

00 

1(5) = J (5 - 5, + i5j) , (1) 

—00 

where g(u) is a functional of the equilibrium velocity distribution 
and , in the initial value problem, represents a Landau or gyro 
resonant normalized velocity. 

The prominence of Maxwellian velocity distributions has led to 
the tabulation^ of the plasma dispersion function Z(?) , defined as 

OO 

=472/ ■*“ (5j > 0) , (2) 

^ 'Lt, 

—00 

and as the analytic continuation of this for ? < 0 . The essential 

2 ^ 

singularities of exp(-u ) at u = ± i «> preclude the computation of (2) 
by contour integration using an infinite semicircle and Jordan's lemma. 
Thus a closed form for the quadrature defining z(5) does not exist, 
and its utilization is awkward. 

It is our purpose here to analyze and compare two simple approxi- 
mations to Z(|) and Z^(|) = dz(§)/d| . The first approximation is 
based on the utilization of resonance velocity distributions, and is 

derived in Section 2. The second approximation is given in Section 3 

2 

and was proposed by Fried et al . Section 4 compares these two approxi- 
mations. Section 5 applies them to Landau and whistler waves, and 
Section 6 discusses the results, commenting on the possible improvement 
of the Fried et al . approximation. 
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2. THE RESONANCE DISTRIBUTION APPROXIMATION 


The Maxwellian velocity distribution 

1 


. 1/2 

(2jt) V 


exp- p 

0 \ 


i) 


( 3 ) 


00 

/ 


<v^>^ = / dv F^(v) 


00 

■/ 


dv fj^(v) = 1 , 


3 th 

is sometimes approximated by the n -order resonance distribution 


■'tot'') = 


[2Vg(2n-3) 

2rt (2n-2); 


|^v^+ (2n-3)vgJ’^ 


/ 2v 2 

\n = ''0 


00 

, / IV F„„(V) = 1 . 


( 4 ) 


In terms of the complex v-plane, the essential singularities of 

"fch 

F (v) at V = ± i <» are simulated by two n^ -order poles at 
^ 1/2 

V = ± i(2n-3) Vq . These poles tend to v = ±i<» as n-*«>^ and , 
as shown in Appendix A, the Maxwellian distribution is retrieved in the 
limiting process: 


llm F^Jv) - F^(v) . 

n -« 00 

This result suggests that z(|) might be approximated by 

CO 





du 

(u^+y^)"(u-?) 


(Ij > 0) 


( 5 ) 


(6) 


wh ere y 

n 


(n- 3 / 2 )^^^ and 


n> 2 


The quadrature can now be performed 


2 



by contour integration. Closing the path with an infinite semicircle 
in the lower complex u-plane, we obtain > O) 



du 


(u +y^)’^(u-?) 


2ni 

(n-l)’. 


, . d”-l 

lim r 

j n-1 

u =-iy du 
n 


1 

(u-iy^)”(u-5) 


( 7 ) 


so that Z„ becomes 
Rn 


V'5) = - 


( n-1 ) (2n-m-2)’. 


(I, > 0) . 

(8) 


Integration by parts readily shows that Z^(|) = dz(|)/d| can be 

2 2 

obtained from (2) by substituting d[exp(-u )]/du for exp(-u ) . 
Approximating the derivative of the Maxwellian by the derivative of the 
resonance distribution amounts to term by term differentiation of (8). 
The resonance distribution approximation of Z^(|) is then 



(?) 


n , vm-1 

(n-lV. -^ (2n-m-l)-. 

(2n-2)’. ^ (n-m): \m+l 




0 ) 


( 9 ) 

We note that it is sufficient to find an approximation to z(|) 

and Z^(§) in one of the four quadrants of the §-plane because of the 

/ N 4 

symmetry properties of Z(§; and its derivatives, 




( 10 ) 


and the knowledge of their analytic continuations into the lower 
5-plane, 


+ 21 (-1)" Hj5)d=q>(-5^) 




5 


(I, < 0) , 

( 11 ) 



where the Hermite polynomials H (|) satisfy the recurrence relation 


n 


“n (5) = 2? - «'.i (5) 


Hq(I) = 1 . 


h 



3. THE TWO-POLE APPROXIMATION 


Instead of approximating the Maxwellian distribution and then 

2 

computing its Hilbert transform, as done in Section 2, Fried et al . 
have approximated Z and h' directly, obtaining for ^ ^ 





1 

a*(a*+§) ' 


a 


1/2 

0.55 + 1 


a + 1 a 
r 1 


( 12 ) 


and 



~ 1 ^ 1 
_b(b-|)^ b*(b*+5)^_ 


i=0.45+i0.86=b +ib 
b r 1 


(15) 


In the lower half-plane, < 0 , Equation (ll) is used to analytically 

continue Z„ and 7i' . 

F F 

Fried et al . have derived this form of Z„(5) by requiring that the 

" r 

two-pole approximation displayed the symmetry properties and asymptotic 

behavior of z(§). The imaginary part of 1/a was obtained by imposing 

the condition z(§=0) = Zj,(§=0) , and the real part of 1/a was chosen, 

with an "eyeballing" procedure, to minimize | Z-Z^,] . In Section 6 we 

shall comment on other possible criteria to choose the value of 1/ a . 

The derivation given in Appendix B, or the substitution in (l) of 

g(u) by the f„(u) given below, shows that the form of the (velocity) 

distribution f„(u) implicitly utilized in the Fried et al. approxi- 
f 

mation of Z(|) is 


i^(a2+a2) . 

r 1 1 

(l .21-Ht)rt^^^ [(u+a )^ + a^][(u-a )^ + a^] 

' r 1 r 1 
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Whereas the essential singularities of the Maxwellian at u = ± i “ 

are simulated by two n^^-order poles at u = ± i(n-5/2)^^^ in the 

resonance distribution method, the approximation given by Fried etal. 

implies^ the use o±^^ velocity distribution with four single poles located 

at u=±a ±ia. (or u=±b ±b. in the case of z') . 

r 1 r 1 F 
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k. COMPARISON OF THE TWO APPROXIMATIONS 
The plasma dispersion function z(§) represents the Hilbert trans- 


form of 


3t 


(15) 


Similarly, the approximations Z (§) and Z (5) are Hilbert trans- 
' ^ Rn 


forms of fp(u) given in (Hi-), and 


(2y ) 


2n-l 


[(n-1)’.]^ 1 

2n-2): / 2 . 2sn 

(u +2yJ 


( 16 ) 


These normalized velocity distributions are depicted in Fig. 1, with 

n = 2 , for u > 0 . We note that f is a very good approximation 

to f,, for small values of u whereas f„ is a better asymptotic 
M Rn 

approximation to fj^ . Because of the type of weighing imposed by the 
numerator of the integrand in ( 2 ) near u « § , we- expect to be 

a better approximation to Z than Z . 

The actual comparison between Z and Z , and their approximations 


( with A = F , R2 and R6 , is given in Figs. 2-4 by 

'a 


Z. and Z^ 

A A - (') (') 

representing the complex error A(5) = Z^ ''(§)- Z' (5) along three 

distinct paths in the first quadrant of the complex §-plane. It is not 

necessary to explore the other quadrants of the |-plane because of the 

symmetry properties outlined at the end of Section 2. 

These figures show that, except for large ||| , when simpler 

asymptotic expressions may be used to simulate the plasma dispersion 

function and its derivative, the approximations suggested by Fried 

et al . are far superior to the approximations based on the resonance 

distribution. As the illustrations suggest, it would be necessary to 

go to high values of n to improve on the simulation suggested by 

Fried et al . However, since the resonance approximation has n terms 

th 

and involves poles of up to the n -order, it is clear that its use 
for large n becomes awkward and, in the light of the Fried et al . 
approximation, unjustified. 
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FIG. 3. Errors A(|) = Z - Z and z' - z' for A = F , R2 and r 6 along | 



5. APPLICATION TO LANDAU AND WHISTLER WAVES 


The dispersion relation of electron Landau waves propagating as 
exp i(oDt-kz) in a Maxwellian plasma is^ 


K 


= z'(§) 


(■■?) 


(17) 


where we have 


K = 


kv . 


CD 


W = 


CD 

oT 


v2 

t m 

e 


03 IS the electron plasma frequency; k is Boltzmann's constant, and T 
P ^ 

and m^ are the electron temperature and mass. 

In Fig. 5 we plot the solutions of this dispersion relation lor 

real frequencies^ K = K(w = corresponding to the lowest order root, 

and compare them with the curves obtained by substituting z' and z' 

r Rid 

for z' . We use the expressions given in ( 9 ) for n = 2 , and ( 15 ) 
when ^ 0 ^ and 


Z^(C) = 5 exp(-l^) (A = RZ,F) ( 18 ) 

for < 0 . (in Fig. 5 we have •) difficulties 

associated with 1§| as W -* 1 , the curves obtained with lor 

W < 1.08 are extrapolations. For reference, the solution of the Bohm 
and Gross dispersion relation, 

= 1 + I , (19) 


IS also depicted. 

The curves show that the two approximations yield roughly equi- 
valent results: z' is more adequate to compute K , whereas Z^ 

F 1 

gives better agreement with . We note that the resonance approxi- 
mation was used for rather large values of thus explaining its 

relative success (see Section 4). For applications involving small 
111 , the advantages of the Fried et al . approximation become clearer, 
as we shall see for the whistler case. 
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The dispersion relation of electron whistlers propagating as 

exp i(tut-kz) along the static magnetic field in an isotropic Maxwellian 

, . 6 

plasma is 


where we have 



w2 


1 + _E_ z(|) 

WKp 



K 


kc 

6D 





( 20 ) 


n is the electron cyclotron frequency defined by the static magnetic 

field; c is the speed of light in vacuum, and to and v were 

P t 

defined in relation to ( 17 ). 

After choosing W and p , we solve the dispersion relation for 

real frequencies, K = k(w = W^) , and plot the results corresponding 

to the least damped root of (20) together with the solutions obtained 

by substituting Z_ and Z for Z . In Figs. 6-8 we have used 

W =5 values of p corresponding to electron temperatures of 1, 

P 

10 and 100 eV. The approximations used for Z when > 0 are Z^, , 
given by (l2), and , defined by (8) for n = 2 . When 0 

have adopted the following analytic continuations 


Z^(§) = exp(-5^) (a = R2,f) . ( 21 ) 

Here we find ^ W < 1 and W > 1 . For compari- 
son, these figures also show the solution of the cold whistler dispersion 
relation, 



(22) 


Because in the region of cyclotron resonance (w ~ l) we have 
1§I ~ 0 , we expect the results obtained with the resonance method to 
be poor approximations of the exact solutions in that domain (see 
Section 1^) . Indeed, the curves show that whereas Z^, yields a good 
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15 


plasma temperature is 
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PIG. 7. Brillouin diagram of whistler waves using (— ) the plasma dispersion function; ( — ) the 
Fried et al . approximation (l2); (•••) the resonance approximation (8) with n = 2 . The 
solution of the cold plasma dispersion relation is denoted by (— •). The electron plasma 
temperature is 10 eV . 






approximation^ the solutions based on the resonance method yield poor 

results that worsen with an increase in the electron temperature. In 
particular^ we find a discontinuity in the results for lim (w = 1 ± 6) 

6-*0 

which is brought about by the combination of the error Z (§)- z(§) at 

R2 

§ = 0 , and the adopted analytic continuation (2l). This discontinuity is 
removed if we utilize (8) both for > 0 and | < 0 , i.e. if we 

analytically coritinue the algebraic expression (8) disregarding the 
characteristics of the function that (8) is trying to simulate. However^ 
this procedure, although removing the discontinuity, shall bring about 
larger errors than the adopted (21 ) as w(> l) increases. 
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6. DISCUSSION 


The results presented in Section 4 , showing the superiority of the 

two-pole approximation of the plasma dispersion function with respect 

to the resonance approximation, are confirmed by the applications given 

in Section 5. It should be noted, however, that the apparent simplicity 

of both approximations, equations ( 8 ) and (l 2 ), is deceiving. When the 
I 

problem under consideration requires the utilization of the plasma 
dispersion function both in the upper and lower half planes of its 
argument > 0 and 0 , as e.g. in the whistler problem for 

W < 1 and W > 1) it is necessary to utilize the analytic continuation 
given by (ll). The algebraic simplicity of the approximations is then 
lost . 

The criterion adopted by Fried et al . to obtain the two-pole 

approximation was described in Section 3. Accepting the general form 

(12) of the approximation, which satisfies the symmetry properties and 

asymptotic behavior of z(§) , we are left with the choice of the complex 

parameter 1 /a = a . Imposition of the condition z(? = O) = Zy(| = O) 

determines a, = . But, if instead of minimizing [z -zj by 

"eyeballing", is chosen to minimize the square relative error 

/ll - Zj/Zl^ d§ over the most critical path^, i.e, the real axis (here 

taken between | =0 and 4 ), we find that the optimum value of a is 

r ^ r 

« 0.69 and not O.55 as proposed by Fried et al . Also, the condition 
z (5 = 0) = Z (? = 0) might be discarded when Z is not utilized for 
small arguments. Then, minimizing the square relative error between 
I =0 and 4 for a = 0.69 > we find that the optimum value of a. 
is ?» 0 . 85 . Figure 9 depicts the relative errors A = (z-z)/|z| for 
Z! = Z_ given by (l 2 ), and Z = Z where Z^ is obtained from (l 2 ) 
by putting a = 1/a = 0.69 + ^ 0.85 . find that Z^q yields smaller 

errors than Z^ for > 0.4 . 
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(Zp-Z)/|Z| 



KkP 
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FIG. 9. Relative errors A(§) = (z^-Z)/|z| for A = F , i.e. (l2) with 1/a = 0.55 + i 
and A = FO, i.e. (l2) with 1/a = 0.69 + i 0.85 , along 5. = 0 . 
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APPENDIX A 


THE MAXWELLIAN DISTRIBUTION AS A LIMIT 


To demonstrate (5)j we use Stirling's formula to write 

t(n-l)'.? . vl/2 

2 n-a) ~ " '“2 


lim 


(A.l) 


and note that 


lim 


|[2vg(2n-3)l^^] 


2n-l^ 


n -* “/2jt[v +(2n-3)v^] 


= lim 

n -» 00 


1/2 1 o 2 n 

QV-v \ 2 vq ^ 


-n 


jrVo(2n) 


. (a. 2) 


Substitution in (4) yields 


lim = lim 

n “ n " 


(2n) 


1/2 


1 + 


2v' 


2 nj 


-n 


= V''> 


(A. 5) 
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APPENDIX B 


DERIVATION OF f„(u) 
F 


To obtain fp(u) we solve the integral equation 


fp(u) 


u = 


a /2 ^ -1 

2 


, ? > 0 . (B.l) 


Denoting the positive frequency part of f„(u) by fp(u) , that is 


4( 


OS 

' 2n J 


ds f_(s) exp isu 

r 


0 




(B.2) 


1 (s) = / du f (u) exp(-isu) , 


“F'“' J F 

.dOS 


we can write 


a 

/ 




Zp(§) = J du = 1 2n fp(§) 


(5^ > 0) . 


(B.5) 


Because f„(s) is the Fourier transform of f„(u) and we expect this 
velocity distribution to be real and even (it simulates the Maxwellian 
distribution) it follows the f_(s) will also be real and even. 
Combining (l2), (B.2) and (b.j) we find 


2 » 4 ( 5 ) = 


1 

a(^-a) 


00 

a*(5+a*)] = j 
0 


ds f_(s) exp IS 5 
F 


( 5 ^ > 0 ) , 

(B.4) 


SO that introducing p = - i§ we obtain 
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1 ,la(p+ia) 


CO 

1 .la*(p-ia*) J 
0 


ds fp(s) exp(-ps) (p^ = > 0) 


(B.5) 


We retrieve fp(s) , for s > 0 ^ by inverting this Laplace 
transform: 


v«) 


exp( -las) exp ia^s 
1.1 a 1.1 a* 


1/2 


(B.6) 


= cos (a s) exp a s + — - sin (a s) exp a.s (s > O) . 

' r ' 1 1.1 ' r 1 ' ' 


The desired velocity distribution is then obtained by inverting the 
Fourier transform^ and recalling that ^p(s) is a real and even function 


of s 




00 


ds f„(s) exp isu 
£ 


00 

0 


ds i_(s) cos (su) 

r 


2jt 


/ \2 2 ( \2 2 

( u+a ; +a . ( u-a ) +a 

' r' 1 ' r 1 


(B.T) 


2.2n 


1/2 


u + a 


u-a 


( 2 / ^2 2 
(u+a^j + (u-a^) + a^ 


1/2 


Noting that 1.1 |a^| = rt a^ , we finally have 


4(a^ + a^) 
fJu) = ^ ^ 


(1.21-Ht)/^^ [(u+a^)2 + a^][(u-a^)2 + af 


with 


a = 


2.2 


r 1.21+it ’ 


i 1.21+rt 


(B.8) 
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